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ABSTRACT 

We have used the Atacama Pathfinder Experiment 12 m telescope (APEX) to carry out an absorption study of submillimeter 
wavelength rotational ground-state lines of H 35 C1, H 37 C1, I3 CH + , and, for the first time, of the SH + radical (sulfoniumylidene or 
sulfanylium). We detected the quartet of ground-state hyperfine structure lines of SH + near 683 GHz with the CHAMP+ array receiver 
against the strong continuum source Sagittarius B2, which is located close to the center of our Galaxy. In addition to absorption from 
various kinematic components of Galactic center gas, we also see absorption at the radial velocities belonging to intervening spiral 
arms. This demonstrates that SH + is a ubiquitous component of the diffuse interstellar medium. We do not find clear evidence for 
other SH + lines we searched for, which is partially due to blending with lines from other molecules. In addition to SH + , we observed 
absorption from H 35 C1, H 37 C1, and l3 CH + . The observed submillimeter absorption is compared in detail with absorption in 3 mm 
transitions of H 13 CO + and 0C3H2 and the CO 7 = 1-0 and 3-2 transitions. 
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1. Introduction 

1.1. (Sub)millimeter spectroscopy of diffuse interstellar 
clouds 

Starting in the late 1930s, studies of molecular absorption toward 
diffuse and translucent clouds have traditionally been in the 
realm of optical and, later, ultraviolet (UV) and infrared (IR) 
spectroscopy. Over the last 15 years, observations at radio and 
(sub)mm wavelengths have greatly added to our knowledge of 
the chemistry o f the diffuse interstellar medium (ISM) (see, e.g., 
the review of |S now & McCall 2006). Measurements toward 
strong continuum sources are a powerful means to detect spectral 
lines that have their lower level in (or near) the rotational 
ground state of a molecule and may be the only way to detect 
a number of species in low density environments or even at 
all. In particular, the large rotational constants of light hydrides, 
resulting in a rotational spectrum starting at submillimeter or far- 
infrared wavelengths, together with substantial electric dipole 
moments, make collisional excitation negligible compared to 
radiative decay. This leaves substantial population only in the 
lowest energy level, which can be split in the case of fine 
structure (fs) or hyperfine structure (hfs). 

Two particularly fruitful avenues of research have emerged; 
First, choosing extragalactic sources with strong mm- and cm- 
wavelength continuum emission as background objects, H. Liszt 
and R. Lucas have studied the abundances of a variety of di- 
, tri-, and even polyatomic molecules toward up to 20 lines 



of sig ht with often surprising results (see, e.g., iLiszt & Lucasl 
2002, and references therein). A second possibility is using as 
a background the thermal continuum from star-forming regions, 
(mostly) free-free emission at cm- and longer mm-wavelengths 
and dust emission at shorter wavelengths. Particularly interesting 
lines of sights are those toward regions M and N in the 
Sagittarius B2 Giant Molecular Cloud (GMC), which allow 
sensitive absorption spectroscopy of the ISM surrounding the 
Galactic center (GC) near which they are located at a distance 
of ~ 8 kpc dReid et al.l 12009) as well as of intervening 
GMCs. For a long time, absorption between -120 km s _1 
and Sgr B2's local velocity range ( > +50 km s" 1 ) has been 
observed in the HI 21 cm line (see iGarwood & Dicke"vl fl989. 
and references therein) as well as in (mostly) mm- wavelength 
lines from quite a number of molecules (se e iLinke et aTl 
I1981L lNvman|[l98l iGreaves & Wimamsl[T994 and references 
therein). Intervening spiral arm absorption toward other bright 
inner Galaxy (sub)mm continuum sources has been found as 
well, most prominently towar d the central radio/(sub)mm source 
Sgr A* and W49N (see, e.g.. ILiszt et al.ll97^;lLinke etail l98 1 ; 
| Coxetal.lll988t IGreave s & Williams 1994; lPlume etalJ 12004: 
Poleham pton et al.ll2005bl) . 



At swbmillimeter wavelengths, numerous molecules have 
been detected in absorption via their ground-state transitions 
(until recently mostly) toward Sgr B2 (for summaries see 
Menten 2004; iLis et all 120091) . In particular, the ho - l i 
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transition of ortho-H^O and 



Hj 8 was detected with the 



Submillimeter Wave Astronomy Satellite dNeufeld et al.l l2000): 
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this H} 6 line, near 557 GHz, is responsible for much of the 
Earth atmosphere's shorter-wavelength submillimeter opacity. 

In fact, atmospheric absorption by the rotational spectrum 
of H2O (and to a lesser extent of O2) makes the ground- 
state lines of many light hydride species unaccessible from 
the ground. Using the Long Wavelength Spectrometer (LWS) 
aboard the Infrared Space Observatory (ISO), which covered 
frequencies > 1.5 TH z (wavelengths < 200 //m), s everal 
hydroxyl isotopologues (Polehampton et al. 2003, 2005a) were 
detecte d, as well as CH2 (Poleham pton et al.ll2005bl) . and para- 
H2D + dCernicharo et alJ 20071) . From the excellent high altitude 
site in the Atacam a desert, recently OH + was found with the 
APEX telescope bv lWvrowski et al.l d2010al) . 

Very recently, a new era in submillimeter absorption 
spectroscopy started with operation of the High Frequency 
Instru ment for the Far In frared (HIFI; see Ide Graauw et alJ 
120 lOt) aboard Herschel dPilbratt et al.l I20101) . Utilizing the 
superb sensitivity of Herschel's 3.5 m aperture combined 
with HIFI's high spectral resolution heterodyne receivers 
operating in a low thermal background space environment, 
revolutionized observations of light hydrides. (Mainly) in the 
framework of the Guaranteed Time Key Projects Herschel 
observations of EXtra-Ordinary Sources (HEXOS) (PI: E. A. 
Bergin) and PRobing InterStellar Molecules with Absorption 
line Spectroscopy (PRISMAS, PI: M. Gerin) observations were 
made of H 2 + , H 3 + , H 2 and its isotopologues, H 2 C1 + , C 3 , 
CH, and CH + , as well as NH and N% while NH + has so 



far remained undetected (Ossenkopf et al] 1201 Ot ISchilke et af 



2010; Gerin et al. 2010; Wvrowski et al. 


2010b; Neufeld et al. 


2010a: iLisetalJ 


2010b; Melnicketal. 


l20ld IComitoetal 


2010; Lis et al. 2010a; Phillips et al. 2010 


;lNeufeld et al.ll2010b; 


Sonnentrucker et al.| [20 lOt JMookeriea et alJ |2010t lOin et alJ 


2010tlFalgarone et al.ll2010HPersson et alJl2010h. 



High quality optical/UV observations, needed, e.g., for 
meaningful abundance determinations, mostly sample relatively 
nearby clouds since they require optically bright background 
stars, typically members of OB associat ions (with HP catalogue 
numb ers) within a few kpc (see, e.g. lGredel|[T997t IStahl et al] 
120081 and references therein). In contrast, being much less 
affected by interstellar extinction, radio/submm and near 
infr ared absorption obser vations, e.g. of the key H^ molecular 
ion dGeballe & Okal 19961) . are capable of probing diffuse clouds 
throughout the Milky Way. 

1.2. Interstellar SH + 

The ionized hydride SH + (sulfoniumylidene or sulfanylium) 
potentially is a probe of comets, photon-dominated regions 
(PDR s), shocked in terstellar gas and diffuse molecular clouds 
dHorani et al.l [l985) . Toward the latter, optical absorption 
measurements have shown for a very long time large abundances 
of another ionized molecule, methylidynium, CH + , the first 
int erstellar molecule, detected b y lDunhaml (fl937i) and identified 
by [Douglas & Herzberg (Il94ll) . CH + 's abundance is orders of 
magnitude higher than predicted by othe rwise successful models 
(see, e.g., Ivan Dishoeck & HacS 1986). This has been a long- 
standing puzzle to astrochemists, who invoked shock wave 
scenarios to over come the formation reaction's endothermicity 
dFedermanl!982l) . 



Milla r et al J d 19861) and also iPineau des Forets et al] d!986l) 



have extended chemical models by including sulfur chemistry 
and predicted efficient SH + formation to occur in shocked diffuse 
clouds and that the abundance of SH + can approach that of 
CH + in moderate velocity magnetic shocks, but is negligible 



for a non-magnetic shock. This, in principle, could make the 
CH + /SH + abundance ratio a good shock indicator. 

Millar & Hobbs] dl988l) made sensitive searches of an optical 
wavelength electronic transition of SH + from the A 3 n~X 3 £(0, 0) 
band at 3363.49 A toward the well-studied diffuse cloud in the 
direction of £ Oph, which shows significant CH + absorption. 
They were only able to put an upper limit on the SH + column 
density of 1 x 10 13 cm 2 , which is a factor of 3 low e r than 
the measured CH + column density. Magnani & Salzerl dl99ll) 
searched for this line and another one at 3366.08 A toward a 
sample of 13 stars with known diffuse molecular clouds in front 
of them and were only able to determine comparable 3<x upper 
limits on the SH + column density between 0.4 and 1 x 10 13 cm -2 
and upper limits on the [SH + /CH + 1 abundance ratio between 
0.12 and 2. An earlier search by iMagnani & Salzerl (l989) 
toward a sample of diffuse and high-latitude clouds resulted in 
a factor of a few coarser limits on the SH + abundance. High 
(Galactic) latitude clouds belong to the so-called translucent 
clouds. 

Translucent clouds, which were recognized as a class 
by Ivan Dishoeck & Blackl d 19861) have physical properties in 
between those of diffuse and dense dark clouds. In particular, 
they have higher visual extinctions, Ay , (1-a few mag) than 
diffuse clouds (< 1 mag) and, consequently, a higher ratio of 
molecular to atomic gas. 

Other astrochemical interest in SH + derives from the fact 
that, given a possible formation route via doubly ionized sulfur 
(S +2 ), it is predi cted to be a dia g nostic of X-ray dom inated 
regions (XDRs) dAbel et al.ll2008l) . IStauber et al] d2007l) report 
an observation of the SH + N = 1 - 0, J = - 1 , F = 0.5 - 1 .5 
line near 346 GHz toward the high-mass young stellar object 
(YSO) AFGL 2591, which is surrounded by an XDR. Since this 
line is blended with a 34 S02 line, they were only able to quote 
an upper limit to the SH + line's intensity. Again very recently, 
Herschel/HIFI observations led to the discovery of emission in 
the 526 GHz I2 — 0i transition of S H + toward the high-mass 
protostar W3 IRS 5 dBenz etal]|2010l) . 

In ^2] of this paper we review the laboratory spectroscopy 
picture for both SH + and 13 CH + . In ^3] we describe our 
submillimeter observations. Our results and the analysis 
procedure used to derive primary line parameters of the lines 
are presented in §|4] In {(5]we discuss the derived line parameters 
in an astrochemical context. 



2. Rest frequencies 

2.1. Accurate SH + frequencies 



I Savage et al.l ([2004) performed the only measurements of SH + 
with microwave accuracy. They obtained rest frequencies for 
four hfs lines from the X 3 2T = 1-0 state, one near 346 
and th ree near 526 GHz. A previous study bv lHovde & Sav kallv 
(1 19871) employed laser magnetic resonan ce to record rotational 
spectra. Recentlv. lBrown & M iiller (2009) used these and further 
rovibrational (infrared) data on energy levels of SH + in and 
between the v = and 1 vibrational states of the X 3 X~ 
electronic ground state to perform a single, weighted least- 
squares fit in order to determine an improved set of molecular 
parameters for this molecule. They used their results to calculate 
the rotational spectrum of SH + in these vibrational levels up to 
the N = 4-3 transition. The predictions presented in the catalog 
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Fig. 1. Detail of the energy level diagram of SH + . The solid 
arrows mark transitions observed in the course of this work. The 
thickness of the arrows indicates the relative strengths of the 
transitions. The hyperfine splitting has been exaggerated slightly 
for clarity reasons. 



sectiorQ o f the Colog ne Database for Molecular Spectroscopy, 
CDMS. (iMtlller et all l200ll 120051) are based on the data from 
ISavage et al.1 ( 2004 ), on weighted av erages of the experimental 
data from lHovde & Savkailvl (Il987|) extrapo lated to zero field 
taken from the fit of iBrown & Mt iller (2009) and all available 
infrared data (up to v = 4 - 3). These predic tions are essentially 
identical to those in iBrown & Mulleil ( 2009). in particular when 
the uncertainties are taken into account. The dipole moment of 
1.285 D in the ground vi brational state has been ta ken from 
an ab initio calculation bv lSenekowitsch et al.l ([T985). The rest 
frequencies from the CDMS catalog are given in TableQ] Fig.Q] 
presents the lower part of the rotational energy level diagram of 
SH + . 



2.2. The J = 1 - transition of 13 CH + 

The X 1 E + 7 = 1-0 rotational ground-state transition of the main 
isotopologue of methylidynium, 12 CH + , near 835.2 GHz, cannot 
be observed from the ground for objects with Galactic radial 
velocities becau se of a strong telluric O? line approximately 
1.1 G Hz lower dGolubiatnikov & Rrup nov 2003j iDrouin et al.l 
120101) . Observations of the 7 = 1 - line of 13 CH + at an ~5 GHz 
lower frequency are however possible. 

Falga rone et al.l d2005l) reported its detection toward the 
bright star-forming region G10. 6-0.4 and evaluated from their 
observation the rest frequency as 830132 + 3 MHz. Two 
alternative values for the rest frequency, 830107 + 1 or 830193 ± 
4 MHz, remained viable due to ambiguities in ascribing 
the observed feature to LSR velocity components possibly 
associated with this source. These were disc arded on the basis 
of the poor agreement with the prediction of iPearson & Drouinl 
(2006) who had reported a transition frequency for 7=1-0 
of 12 CH + and predicted the corresponding 13 CH + rest frequency 
as 830134.3 MHz. This value was based on a combined analysis 
of their frequency and of optical CH + data of four isotopologues 
employing an isotopically scaled Dunham analysis. 



Our observations of 13 CH + , however (see indicate a rest 
frequency of around 830.2 GHz if we assume that this cation has 
an absorption profile similar to those of SH + , HC1, H 13 CO + , or 
C- C3H2. 

Miiller (2 010|) performed an analy sis similar to the one 
published by IPearson & Drouinl d2006l) in order to evaluate 
rotational data for isotopologues of CH + independently. 
Applying appropriate parameters which describe the breakdown 
of the Born-Oppenheimer approximation, he obtained a rest 
frequency for 13 CH + 7=1 - com patible with the present 
observations. The reviewer of iMullerl (I2010h pointed out new 
laboratory measurements of the 7 = 1-0 transitions of 
12 CH + , 13 CH + , an d 12 CD + which were submitted subsequently 
bv lAmanol d2010h . Not only was the value for 13 CH + in quite 
good agreement with the estimate from our observations, but it 
was also shown that the previously reported value for 12 CH + 
was i ncorrect by sev eral tens of mega hertz. I ncludi ng the data 
from lAmanol d2010l) in his analysis, Mullet] d2010l) calculated 
transition frequencies for several isotopologues of CH + which 
now should be quite reliable and which are in the CDMS catalog. 
13 CH + 7=1-0 value is again given in TableQ] 

The dipole moment of 1.683 D for CH + in the ground 
vibra tional state has be en taken from an ab initio calculation 
by ICheng et al.l d2007l) : the change in dipole moment upon 
substitution of 12 C with 13 C is assum ed to be small. A value o f 
1.7 D was used in the CDMS catalog dMiiller et al.ll200lll2005l) . 

2.3. HCI, C-C3H2, H"CO + , and CO 

The rest frequencies for the 7=1-0 transitions of both 
H 35 C1 an d H 37 C1, which ha ve hfs, were taken from the JPL 
catalog0( Pickett et al. 1998). These frequencies were published 
by iDe Lucia et all d 197 lb . The dipole moment of H 35 C1 for 
v = 0, 7 = 1 has been determined as 1.1086 (3) D by 
Ide Leluw & Dvmanusl dl973l) : the dipole moment of H 37 C1 is 
expected to be only marginally different. See Table Q] for the 
frequencies of HCI and the other species considered. 

The estimated accuracy for each of the three hyperfine 
components is 100 kHz, sufficient for the present analysis. 
It should be mentioned in this con text that improved rest 
frequencies have be en published by iKlaus et a D d!998l) and 
ICazzoli & Puzzarinil d2004l) employing Lamb-dip spectroscopy. 

The C-C3H2 frequency data were also taken from the JPL 
catalo g. The 2^ - lot transition was measured by Vrtile k et al.l 
dl987h : the reported uncertainty is smaller than the value used in 
the JPL catalog. 

The H 13 CO + and CO rest frequencie s were taken from the 
CDMS catalog dMuller et al.l l200 1 . 2005). The rest frequency of 
the 7=1- transition of H 13 CO + is f rom an astronomical 
observation by Schmid-Burg k et al.l d2004l) . both CO transitions 
were from Lamb-dip measurements by IWinnewisser et al] 
(fl997h . 



3. Observations and data reduction 

3.1. Overview 

We have used the Atacama Pathfinder Experiment (APEX) 12 m 
telescope for observations of several sources in a number of 
submillimeter wavelength transitions of SH + whose accurate 
frequencies have recently been determined (see these 



1 http://www.astro.uni-koeln.de/cdms/entries/ 
http://www.astro.uni-koeln.de/cdms/catalog/ 



also 



- website: http://spec.jpl.nasa.gov/ftp/pub/catalog/catdir.html 
also http://spec.jpl.nasa.gov/ 
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Table 1. Spectroscopic data and beam-averaged continuum brightness temperature assumed for the LTE modeling of the absorption 
lines toward Sgr B2(M). 
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Position 6 


1 ' 

frequency 


S 


Ei 


F eff 


S eff 


^cont 


reference 7 








(",") 


(MHz) 




(K) 


(K) 




(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


C-C3H2 2i2— lot 






2.6,-2.0 


85338.893(14) 


1.50 


0.0 g 


0.95 


0.78 


5.9 


JPL,/i 


H 13 CO + 1-0 






2.6-2.0 


86754.2884(46) 


1.0 


0.0 


0.95 


0.78 


5.9 


CDMS, i 


CO 1-0 






2.6,-2.0 


115271.2018(5) 


1.0 


0.0 


0.95 


0.74 


5.9 


CDMS, j 


CO 3-2 






2.6,-2.0 


345795.9899(5) 


3.0 


16.6 


0.95 


0.72 


11.0 


CDMS, / 


SH+ lo-O, 


0.5 


-0.5 


0,0 


345843.591(128) 


0.2289 


0.0 


0.95 


0.72 


11.0 


CDMS, - 




0.5 


- 1.5 




345929.810(75) 


0.5780 


0.0 








CDMS, k 


H 37 C1 1-0 


1.5 


- 1.5 


0.0 


624964.374(100) 


1.3333 


0.0 


0.95 


0.41 


33.8 


JPL, / 




2.5 


- 1.5 




624977.821(100) 


2.0 


0.0 








JPL, Z 




0.5 


- 1.5 




624988.334(100) 


0.6667 


0.0 








JPL,/ 


H 35 C1 1-0 


1.5 


- 1.5 


0.0 


625901.603(100) 


1.3333 


0.0 


0.95 


0.41 


33.8 


JPL,/ 




2.5 


- 1.5 




625918.756(100) 


2.0 


0.0 








JPL, / 




0.5 


- 1.5 




625932.007(100) 


0.6667 


0.0 








JPL, / 


SH + li-0. 


1.5 


-0.5 


0.0 


683336.08(51) 


0.1890 


0.0 


0.95 


0.41 


30.0 


CDMS, m 




0.5 


-0.5 




683361.98(48) 


0.3778 


0.0 








CDMS, - 




1.5 


- 1.5 




683422.31(50) 


0.9441 


0.0 








CDMS, m 




0.5 


- 1.5 




683448.20(50) 


0.1887 


0.0 








CDMS, m 


13 CH + 1-0 






0.0 


830216.096(22) 


1.0 


0.0 


0.95 


0.35 


47.7 


CDMS, n 


SH + 2,-1, 


0.5 


-0.5 


0.0 


893047.60(95) 


0.2889 


32.8 


0.95 


0.35 


43.5 


CDMS, - 




0.5 


- 1.5 




893073.49(96) 


0.1445 


32.8 








CDMS, - 




1.5 


-0.5 




893107.92(105) 


0.1444 


32.8 








CDMS, - 




1.5 


- 1.5 




893133.81(94) 


0.7225 


32.8 








CDMS, - 



Notes: "Hyperfine quantum numbers as far as appropriate. ^Equatorial offset (J2000) with respect to Sgr B2(M). (0,0) corresponds to (a, <5)j2ooo - 
17 h 47 m 20 s 2, -28°23'05". Torward efficiency. f/ Main-beam efficiency. '"The beam-averaged continuum brightness temperature was derived from 
the depth of the saturated absorption lines in the 3 mm wavelength range, and from a direct measurement of the baseline level at higher 
frequency. 'Catalog used for rest frequency and reference for line measurem e nt, see also section [2] s L o wer level is ortho g r ound state; it is 
2.35 K a bove the zero level IVrtilek et al.l d!987l). ISchmid-Burgk et al l d2004h . - IWinnewisser et all d!997h . ISavage et alj d2004i) . iDe Lucia et all 
dl971h . " iHovde & Savkallvl dl987l) : lBrown & Mulled d2009h . ' lAmanol d2010l) 



included the quartet of ground-state hfs lines near 683 GHz. 
Our observations, described in 33.21 led to the detection of 
these lines in absorption against Sgr B2(M) and (N). In addition 
to absorption from gas associated with the GC, we also find 
absorption at the radial velocities belonging to intervening spiral 
arms, which demonstrates that this species is widespread in the 
Galactic ISM. 

For diverse reasons, our results for other SH + lines were 
not as clear as those obtained for the 683 GHz line, namely the 
N - 1-0, 7 = 0-1 lines near 346 GHz, which are completely 
blended with the CO 7 = 3-2 line and our observations of 
the N = 2 — 1,7 = 1-1 lines were inconclusive. We made, 
however, an unambiguous detection of absorption in the 7=1-0 
transition of the 13 CH + radical near 830 GHz, again from gas in 
the Sgr B2 region and from spiral arm clouds. Absorption in 
that line has previously been reported bv lFalgarone et alj £2005) 
toward G10. 6+0.4, leading to an estimate of its rest frequency 
that is discussed in 32.21 

Finally, we took a wide bandwidth spectrum that covers 
the 7 = 1-0 transitions of both the H 35 C1 and H 37 C1 
isotopologu es of hydrochlor i c acid . The former had been 
observed bv lZmuidzinas et all d 19951) toward a position midway 
between Sgr B2(M) and (N). Since their spectrum only covered 
the LSR velocity range between -10 and +140 km s it does 
not contain information for the spiral arm feature velocities. 
While we confirm the Sgr B2 velocity absorption, no clear 
absorption is found toward the latter in our more sensitive and 
wider bandwidth spectra. In 34.1.41 our upper limits for these 
diffuse lines of sight are compared with the tentative ultraviolet 



detect ion of HC1 absorption toward £ Oph by iFederman et al.l 
(11995b . 

3.2. APEX observations and data reduction 

Our SH + observations were made on 2008 September 19 
and 25 under good to very good weather conditions with 
the 12-m Atacama Pathfinder Experiment telescope, APEX0 
( Glist en et al.1 120061) . The precipitable water vapor (PWV) 
content was between 0.3 and 0.7 mm throughout the 
observations. 

The 683 and 893 GHz lines were measured in the lower 
sideband (LSB) with th e MPIfR-built CHAMP+ receiver array 
( Kasemann et al. 2006). CHAMP+ consists of two modules of 
7 pixels each with one central pixel and the others forming 
a hexagon around it. The one module covers the 350 jum 
atmospheric window and the other the 450 //m window. 
Calibration was obtained using the chopper wheel technique, 
considering the different atmospheric opacities in the signal and 
image sidebands of the employed double sideband receivers. 
The image sideband was rejected to the 10 dB level. During 
the observations the single sideband system temperatures at 
683 GHz were between 1700 and 2600 K. To ensure flat 
spectral baselines, the wobbling secondary was chopped with 
a frequency of 1.5 Hz and a throw of 240" about the cross 

3 This publication is based on data acquired with the Atacama 
Pathfinder Experiment (APEX). APEX is a collaboration between 
the Max-Planck-Institut fur Radioastronomie, the European Southern 
Observatory, and the Onsala Space Observatory. 
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elevation axis. The wobbler was operated in symmetric mode, 
which means that source and off position are interchanged 
between subsequent subscans, which cancels any asymmetries in 
the optical paths. Such observations deliver a reliable estimate of 
the continuum level. We found, however, that the quality of the 
baseline deteriorated with increasing continuum flux density. 

The radiation was analyzed with a new incarn ation of the 
MPIf R Fast Fourier Transform spectrometer (FFTS) dKlein et all 
2006), which provided two 1.5 GHz wide modules set to 1024 
frequency channels each. The resultant channel spacing was 
1.46 MHz, corresponding to 0.49 and 0.70 km s -1 at the 
highest and lowest frequency we observed (893 and 625 GHz, 
respectively). The FFTS modules were operated in series with 
an overlap of 300 MHz to provide a total coverage of 2.4 GHz 
bandwidth. To check the telescope pointing, drift scans were 
made across the continuum of Sgr B2(N). Pointing corrections 
were derived from these measurements. The pointing was found 
to be accurate to within * 3", acceptable given the FWHM beam 
size, 6 B , which is 9" FWHM at 683 GHz. 

The SH + line near 346 GHz was observed with the 345 GHz 
atmospheric window module of the Swedish heterodyne facility 
instrument for the APEX telescope (Vassilev et al. 2008) with 
system temperatures between 320 and 370 K. Two facility 
FFTS backends with 1 GHz bandwidth and 8192 channels 
each were connected to the receiver with a slight overlap to 
reach a total bandwidth of 1.8 GHz. The beam width at the 
observing frequency is about 18". Pointing was again checked 
with continuum scans on Sgr B2(N). 

Observations of the 13 CH + 7=1-0 lines near 830 GHz were 
performed using the CHAMP+ array on 2009, August 1 1 under 
extremely good weather conditions with a PWV content of 0.14 
to 0.2 mm. Two 1.5 GHz FFTS modules were used with overlap 
to reach a total bandwidth of 2.4 GHz. Prior to the integration 
on Sgr B2, the telescope was pointed on the continuum of 
Sgr B2(N) and G10.62-0.38. For the line observations again 
the wobbler was used with a throw of 240" and a chopping 
frequency of 1.5 Hz. 

Finally, observations of the H 35 C1 and H 37 C1 7=1-0 
transitions were made with CHAMP+'s lower frequency module 
in 2009 June under superb weather conditions. 



4. Results and data analysis 

During our three observing campaigns, we have gathered 
data that allow unambiguous identification of the SH + and 
13 CH + molecules toward Sgr B2. In addition, we have made 
observations of the 7=1-0 ground-state transition of both 
the H 35 C1 and H 37 C1 isotopologues of hydrogen chloride. 

In absorption spectroscopy, the column density in a 
transition's lower energy level, N\, can be calculated from its 
observed optical depth, r and line width Ay if its excitation 
temperature, T ex , is known: 



Ni 



±JLU 



tAv, 



(1) 



where A„i the Einstein A coefficient; h and k are the Planck 
and Boltzmann constants, respectively, and v = E u — Ei the 
transition's rest frequency; E u and E\ are its upper and lower 
energy levels, respectively. 

The total column density, N tot , is given by 



where Q is the partition function for the rotation temperature, 

Trot- 
In the case of a line with a single absorption component, 
the optical depth and line width are usually determined by a 
Gaussian fit to the line profile and t is calculated from the line 
to continuum intensity ratio: 



T = -ln(l-|7 L |/r c ) 



(3) 



where 7l and 7c are the observed line and continuum brightness 
temperature, respectively. 

Things are generally not so simple for the immensely line- 
rich Sgr B2 sources, for which practically the whole (sub)mm 
range contains a large number of lines from different species 
in each chosen frequency interval, many of them blended. To 
make things even more complex, many of the simpler species 
show multiple absorption components originating from various 
spiral arms along the line of sight that are blended with each 
other and with spectral components from lines of other species 
at GC velocities. Lines showing hfs, such as SH + are yet more 
diffic ult to analyze. To d eal with this, we employ, as described 
bv lBellocheetal] (l2008h . the XCLASS0 program created and 
developed by Peter Schilke. XCLASS is an extension of the 
CLASS spectral line data reduction program that is part of the 
GILDAS0 software package developed by the Institute for Radio 
Astronomy at Millimeter Wavelengths (IRAM). Using XCLASS 
one may fit all lines in a spectrum simultaneously as described 
in ED 

For each species, XCLASS assumes that the level 
populations are described by a single excitation temperature, 
which we refer to as the rotation temperature, T mt . Different 
species can have different rotation temperatures. Here we divide 
between two cases: First, the critical densities of all the lines 
in the frequency range of interest are lower than the density in 
the modeled emission region. Then we are dealing with local 
thermodynamic equilibrium (LTE), which applies, e.g., for the 
extremely line rich hot cores associated with Sgr B2(M) and (N). 
In particular in the latter, the complex molecules NH2CH2CN, 
C2H5OCHO, and C3H7CN along with the I3 C isotopologues 
of vinyl cyanide have recently been identi fied using XCLASS 
with model temperatures of * 100-150 K (iBelloche et al.ll2008l 
2009; iMiiller et al.l l2008). In the second, very low density case, 
which applies for the absorption lines discussed in this paper, 
7 rot equals the temperature of the cosmic microwave background 
radiation, 7cmb, i-e. 2.728 K. While this might be a good 
assumption for clouds far from the strong Sgr B2 continuum 
sources, e.g., in the spiral arms, assuming 7 rot = 7cmb will 
underestimate column densities for molecules in the vicinity of 
Sgr B2(M) and (N). 

4.1. Modeling of the absorption lines toward Sgr B2(M) 

4.1 .1 . The model for the spiral arm and galactic center LSR 
velocity features 

The SH + (li-0i) spectra obtained with the CHAMP+ array 
consist of the superposition of a number of absorption 
components associated with the dense core Sgr B2(M) and 
several molecular clouds distributed alon g the line of sight 
with d ifferent systemic velocities (see, e.g jGreaves & Williams! 
1994). 



Nw = ^e E '' kT -Q(T mt ), 
gi 



(2) 



4 http://www.astro.uni-koeln.de/projects/schilke/XCLASS 

5 http://www.iram.fr/IRAMFR/GILDAS/ 
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Table 2. Velocity ranges and 12 C/ 13 C isotopic ratios of the 
Galactic spiral arms along the line of sight of Sgr B2(M). 



Velocity range 


12 C/ 13 C 


Location 


(1) 


(2) 


(3) 


-62 


20 


Galactic Center, Sgr B2 


39-25 


40 


Scutum/Crux arm 


22-12 


60 


Sagittarius arm 


8 --9 


20 


Galactic Center 


-13 --50 


40 


3 kpc expanding ring 


< -53 


20 


Galactic Center 



The association of the absorption components' LSR 
velocities with Galactic spiral arms or kinematic features 
associated with the GC located along the line of sight 
is presented in Table [2] It is essentially based on the 
velocity interpretation of H?C O and CS absorp ti on b y 
IWhiteoak & Gardnerl d 19791) and iGreaves & Williams! (|1994l) . 
respec t ively, as well as the analyses of ISofud (|2006) and 
Vallee (2008). Although in the GC direction absorption close 
to Vi sr = could in principle come from any location on the 
line of sight, we favo r a phy sical association with the GC since 
iGardner & Whiteoakl (1 19811) found low values of the I2 C/ 13 C 
ratio (20-30) for this velocity range, consistent with a Galactic 
center origin. 

The blending of these absorption components is further 
complicated by the hyperfine structure of the SH + transition. 
To allow meaningful fitting, we used information derived from 
the absorption spectra of other molecules obtained as part of a 
complete line survey done with the IRAM 30 m telescope in the 
3 mm atmospheric window toward both hot cores Sgr B2(M) 
and Sgr B2(N) fsee lBelloche et al.ll2008ll2009l) . 

The line survey was analysed with the XCLASS software 
(see g4]0 . Each molecule identified in the survey was modeled 
with the following free parameters: source size, temperature, 
column density, linewidth, and velocity offset with respect to the 
systemic velocity of the source, 62 km s for Sgr B2(M). When 
several components were needed to reproduce the observed 
spectrum of a molecule, we fitted a set of parameters for each 
component. The combination of the predicted spectra of all 
identified molecules was done a ssuming that the emi ssion adds 
up linearly (for more details, see lBelloche et al . 200§. 

To model the absorption lines observed toward Sgr B2(M), 
we assumed a beam-averaged continuum brightness temperature 
Tcont of 5.9 K in the 3 mm window, as derived from the 
saturated absorption spectra of HCN and CN (see Table [TJ. The 
modeling of the absorption components was done in the T I0t 
= Tqmb approximation with a flag indicating that the absorption 
spectrum has to be computed against a background emission 
consisting of both the continuum emission and the contribution 
of the other molecules seen in emission, i.e. for each absorption 
component with the equation 



T v = (J v (T ex ) - J v (T mnt )) (1 - e- T >) + C e- T " 



(4) 



with T ex the excitation temperature, T cont the background 
brightness temperature, T"" the emission line spectrum, and t y 
the opacity of the line at frequency v, for a beam filling factor of 
1. 
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6 Where noted, for the components associated with Sgr B2 (around 
+62 km s~'), a r iot > Tqmb was assumed. 



Fig. 2. From bottom to top, panels a to c show absorption spectra 
taken with the APEX telescope with the CHAMP+ receiver 
array. Spectra a and b were taken with the central pixel of 
CHAMP+ toward Sgr B2(M) (see Fig. |3}. Spectrum b was 
obtained toward a position (-20.4", 1.4") offset from Sgr B2(M) 
with the CHAMP+ pixel 2. a: H 35 C1 1-0 absorption spectrum. 
b: SH + 1 1 — 0i absorption spectrum, c: 13 CH + 1-0 absorption 
spectrum. In panels a- c, the red curve is the synthetic spectrum 
described in Sect. 14.11 and the horizontal dotted line shows the 
continuum level, d and e: same as b for the 0C3H2 2iz~ loi an d 
H 13 CO + 1-0 transitions observed toward Sgr B2(M) with the 
IRAM 30 m telescope, respectively. Panel / shows the optical 
depth of the 21 cm HI line vs. LSR velocity as derived by 
the procedure described in §5.1.21 In panels d and e, the red 
curve is the synthetic spectrum of the molecule, the green curve 
is the predicted spectrum including in addition all molecules 
identified in our complete 3 mm line survey so far, and the 
horizontal dotted line indicates the zero level (after baseline 
removal because the observations were performed in position- 
switching mode with a far OFF position, which yields a very 
uncertain baseline level). In all panels, the spectra are plotted in 
main-beam brightness temperature scale, the vertical dotted line 
marks the systemic velocity of the source (62 km s~'). 

That table also lists the 12 C/ 13 C isotopologic ratios we 
assumed. The 12 C/ 13 C isoto pic ratio of the G alactic Center 
components was fixed to 20 (Mti ller et alJl2008h . This value is 
compa tible with the Galactic gradient derived by iMilam et al.l 
(2005, equation 5). This equation was additionally employed 
to fix the ratio to 40 for the Scutum and Norma arm 
components, and 60 for the Sagitta rius arm components, with the 
galactocentric distanc es of [Va llee (2008). The ls O/ 18 isotopic 
ratio was taken from [Wilson & Rood! (1 19941) : 250 for the GC 
components, 327 for the Scutum and Norma arm components, 
and 560 for the Sagittarius arm components. 
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-100 -50 50 -100 -50 50 -100 -50 50 

Velocity (km s 1 ) Velocity (km s 1 ) Velocity (km s l ) 



Fig. 3. a to d Column density profiles computed by adding the Gaussian velocity absorption components listed in Table [3] The 
components with upper limit were not included, e to g Abundance profiles relative to H2. The H2 column density was computed 
from the HCO + column density assuming an HCO + abundance of 5 x 10~ 9 . h to j Abundance profiles relative to H 1. The H 1 column 
density was computed from comparing an H 1 absorption with an H 1 emission profile (see text). 



arbitrarily fixed to 2.7 K. The parameters of the fit to H 13 CO + 
are listed in Table [4] The model predicts too strong absorption 
for the H 13 CO + (1-0) components with v < -50 km s _1 , 
while the corresponding HCO + (1-0) components are well 
fitted (not shown here). The disagreement must result either 
from an underestimate of the 12 C/ 13 C isotopic ratio or from 
contamination by emission lines that are not yet included in 
our complete model. It could in principle also result from larger 
excitation temperatures for HCO + compared to H 13 CO + . 

The fit to 0C3H2 was done with velocity and linewidth 
parameters slightly different from HCO + , although we tried 
to stay as close as possible to the components identified in 
HCO + . The parameters of the fit to 0C3H2 are listed in Table[3] 
The model consists of one component in emission (we did 
not attempt to model the possible wing emission in this case), 
one absorption component to mimic the self-absorption, and 
30 absorption components associated with the diffuse clouds, 
with an excitation temperature of 2.7 K and linewidths varying 
between 3 and 6 km s _1 . 

4.1.2. Modeling the SH + absorption 

Since we considered 0C3H2 a priori to be a better diffuse gas 
tracer than 13 HCO + (but see © and because of the added 
uncertainty of assigning a 12 C/ 13 C ratio to each absorption 
components, we decided to model the SH + absorption spectrum 
based on the parameters derived for 0C3H2, i.e. we used 
the same number of absorption components with the same 
velocities and line widths, except for the component physically 
associated with Sgr B2(M) itself. Again, we fixed the excitation 
temperature of each component to 2.7 K. Since the observations 
were done in wobbler mode, we trusted the baseline level 
and derived the level of continuum emission directly from the 




-100 -50 

Velocity (km s 1 ) 



Fig. 4. Ratio of the column density profiles of CH + and SH + 
shown in Fig. [3] 



We present, in Fig. [2] d and e, the results of our modeling 
for c-C 3 H 2 (2i2-loi) and H 13 CO + (1-0), respectively. The model 
for H 13 CO + (1-0) was obtained by fitting HCO + , H 13 CO + , and 
HC ls O + simultaneously. The fit to the observed H 13 CO + (l- 
0) spectrum consists of four components in emission (one for 
the main emission peak and three to reproduce the wings), 
one absorption component to mimic the self-absorption in 
Sgr B2(Mfl and 34 absorption components associated with 
the diffuse clouds along the line of sight, with linewidths 
varying between 3 and 8 km s _1 and an excitation temperature 



7 Although the fits to the self-absorption features in Fig. [2] d and 
e look relatively good, we emphasize that a uniform r rot = Tqmb 
model is not complex enough to provide reliable physical parameters 
to describe the self-absorption. The parameters listed for the self- 
absorption components in the following tables should therefore be 
viewed with caution 
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Fig. 5. SH + absorption spectra obtained toward Sgr B2(M) with the 7-pixel array CHAMP+ at the APEX telescope. The equatorial 
offset position of each pixel relative to Sgr B2(M) is given in the bottom left corner of each panel. In each panel, the dotted line 
shows the continuum level as measured over the emission/absorption-free channels, and the red curve is the spectrum predicted by 
the r rot = Tcmb model described in Sect. 14. II All spectra are plotted in main-beam brightness temperature scale. 



observed spectra shown in Fig. [5] To compute the synthetic 
spectra, we assumed that the physical parameters of the 
absorbing components are uniform across the CHAMP+ field 
of view and therefore we used the same set of parameters for 
all pixels. The resulting synthetic spectra are overlaid on the 
observed spectra in Fig. [5] The agreement is quite good in 
general, which gives us some more confidence in the way we 
estimated the continuum level. However, significant deviations 
are seen for pixels 1, 3, and 4. As mentioned in §|3] the 
baseline quality of the CHAMP+ array drops with increasing 
level of continuum emission and line contamination from the 
hot molecular core associated with Sgr B2(M). Pixel 4 is the 
most affected, and we believe that the poor quality of the fit 
is mainly due to this problem. As far as the northern pixels 1 
and 3 are concerned, the main absorption component is much 
stronger than the one predicted by our simple model. This may 
be due to spatial variations of the SH + column density across 
the Sgr B2 molecular cloud since the cloud is far from being 
uniform (in particular, the hot, dense core Sgr B2(N) is located 
~ 50" to the north of Sgr B2(M)). If the variations from pixel to 
pixel can be trusted, then we can estimate the variations of the 
SH + column densities in the spiral arm diffuse clouds over the 
footprint of CHAMP + by modeling each pixel separately. This 
exercise yields column density variations typically on the order 
of a factor 2, but in some cases up to a factor 5. However, given 
the low signal-to-noise ratio, the unknown baseline quality, and 
the line blending, these numbers are very uncertain and it is 
unlikely that they really reflect the spatial variations of the SH + 
column density. 

The parameters used to fit the SH + 1 1 — 0i absorption 
spectrum are listed in Tableland the fit to the spectrum obtained 



with pixel 2 is plotted again in Fig. [2J3 for a direct comparison 
with HC1, 13 CH + , c-C 3 H 2 , and H 13 CO + . Although the decision 
to fix the linewidths and velocities to those derived from 0C3H2 
helps a lot to constrain the SH + model, the hyperfine structure of 
SH + and the limited signal-to-noise ratio of the observed spectra 
do not allow to find a unique solution. Using the somewhat 
different set of velocity and linewidth parameters derived from 
HCO + , we were able to find a synthetic spectrum for SH + as 
good as the SH + model derived from 0C3H2 (the parameters of 
this alternative model are listed in Table 0}. Therefore, the SH + 
model described in Table [3] should be considered as a possible 
model rather than the best-fit model. Nevertheless, if we assume 
that C-C3H2 and SH + trace the same diffuse clouds, then it makes 
sense to compare the column densities listed in Table [3] Column 
10 lists the column density ratio SHV0C3H2. It varies by nearly 
two orders of magnitude, from < 0.6 to ~ 33. As a caveat, 
we remind that the SH + model assumes no contamination from 
emission lines of other molecules. If there are contaminating 
emission lines, then the column densities derived for the SH + 
components contaminated by these emission lines are lower 
limits. 

Unfortunately, the SH + 1 — 1 line at 345.9 GHz is heavily 
blended with the spiral arm components of the CO 3-2 line 
seen in absorption (see Fig. |6^). The T mt - 7cmb model 
derived for the 1 1 — 0i line predicts I0-O1 spiral arm components 
in absorption that are relatively consistent with the observed 
spectrum dominated by the CO 3-2 features, except around 
62 km s _1 , the systemic velocity of Sgr B2(M). At 893 GHz, 
there is no clear evidence for the SH + 2i — 1 1 line (see Fig. 0. 
Our T Iot - Tqmb model predicts no significant absorption for 
this higher energy transition (see red spectrum in Fig. |7). On 
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Table 3. Parameters of the best-fit LTE models of SH + , 13 CH + , and c-C 3 H 2 , as well as the alternative LTE model of H 13 CO + fitted 
with the same velocity components as C-C3H2. 











SH + 




13 CH + 




c-C 3 H 2 




H 13 CO 


+ 


n% 


. 1 h 

v hI b 


AV C 


\2q<1 
— 


N e 


N e 


Rl3 f 


R l2 s 


N e 


R" 


N" 






(kms-') 


(km s~' ) 




(cm- 2 ) 


(cm- 2 ) 






(cm- 2 ) 




(cm- 2 ) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


1 


- 


- 




- 


— 


- 


- 


2.4 x 10 14 


- 


- 


- 


2 


64.0 


19.0 




2.0 x 10 14 


1.8 x 10 13 


10.8 


0.53 


1.9 x 10 14 


1.0 


1.4 x 10 14 


0.07 


3 


38.5 


5.0 


■ 


< 2.8 x 10 12 


3.8 x 10" 


< 7.4 


< 0.18 


4.2 x 10 12 


< 0.7 


2.4 x 10" 


< 0.29 


4 


32.0 


5.0 


40 


< 2.8 x 10 12 


< 1.6 x 10" 


— 


— 


4.7 x 10 12 


< 0.6 


7.8 x 10" 


< 0.09 


5 


22.0 


5.0 


60 


1.1 x 10 13 


< 1.6 x 10" 


> 70.6 


> 1.15 


4.1 x 10 12 


2.8 


2.5 x 10 12 


0.07 


6 


17.0 


5.0 


60 


2.1 x 10 13 


9.6 x 10" 


22.1 


0.36 


5.9 x 10 12 


3.6 


2.0 x 10 12 


0.17 


7 


12.0 


6.0 


60 


2.8 x 10 13 


1.4 x 10 12 


19.6 


0.32 


2.7 x 10 12 


10.5 


1.8 x 10 12 


0.25 


8 


8.2 


5.0 


20 


2.5 x 10 13 


2.1 x 10 12 


12.2 


0.60 


8.1 x 10 12 


3.2 


6.0 x 10" 


2.06 


9 


2.5 


5.0 


20 


1.8 x 10 13 


3.5 x 10 12 


5.2 


0.26 


9.7 x 10 12 


1.9 


3.0 x 10 12 


0.30 


10 


-2.5 


4.5 


20 


1.1 x 10 13 


1.9 x 10 12 


5.9 


0.29 


5.4 x 10 12 


2.1 


1.2 x 10 12 


0.46 


11 


-8.7 


5.0 


20 


1.6 x 10 13 


1.6 x 10 12 


9.7 


0.48 


1.4 x 10 12 


11.1 


3.2 x 10" 


2.38 


12 


-15.5 


6.0 


40 


8.4 x 10 12 


2.9 x 10 12 


2.9 


0.07 


3.2 x 10 12 


2.6 


3.6 x 10" 


0.57 


13 


-20.0 


5.0 


40 


8.4 x 10 12 


1.6 x 10 12 


5.3 


0.13 


2.7 x 10 12 


3.2 


5.4 x 10" 


0.38 


14 


-25.0 


6.0 


40 


8.4 x 10 12 


2.5 x 10 12 


3.3 


0.08 


5.6 x 10 12 


1.5 


5.4 x 10" 


0.38 


15 


-29.5 


4.0 


40 


8.4 x 10 12 


8.0 x 10" 


10.6 


0.26 


2.1 x 10 12 


3.9 


3.3 x 10" 


0.63 


16 


-35.5 


4.0 


40 


8.4 x 10 12 


4.8 x 10" 


17.7 


0.43 


2.7 x 10" 


31.5 


1.5 x 10" 


1.38 


17 


-41.3 


3.3 


40 


7.0 x 10 12 


7.3 x 10" 


9.6 


0.24 


4.5 x 10 12 


1.6 


1.6 x 10 12 


0.11 


18 


-45.7 


3.0 


40 


4.2 x 10 12 


8.0 x 10" 


5.3 


0.13 


2.6 x 10 12 


1.6 


5.5 x 10" 


0.19 


19 


-49.0 


3.0 


40 


4.2 x 10 12 


7.7 x 10" 


5.5 


0.14 


3.0 x 10 12 


1.4 


2.9 x 10" 


0.36 


20 


-53.5 


4.0 


20 


9.3 x 10 12 


9.2 x 10" 


10.0 


0.49 


1.2 x 10 12 


7.9 


2.0 x 10" 


2.29 


21 


-58.7 


4.0 


20 


< 2.8 x 10 12 


7.0 x 10" 


< 4.0 


< 0.20 


3.2 x 10 12 


< 0.9 


6.6 x 10" 


< 0.21 


22 


-69.8 


6.0 


20 


9.9 x 10 12 


2.1 x 10 12 


4.8 


0.23 


1.6 x 10 12 


6.1 


5.1 x 10" 


0.96 


23 


-76.0 


4.5 


20 


9.0 x 10 12 


1.6 x 10 12 


5.7 


0.28 


3.3 x 10 12 


2.7 


5.1 x 10" 


0.86 


24 


-80.1 


4.0 


20 


7.0 x 10 12 


1.1 x 10 12 


6.3 


0.31 


4.8 x 10" 


14.6 


2.4 x 10" 


1.43 


25 


-84.3 


5.0 


20 


1.3 x 10 13 


2.5 x 10 12 


5.0 


0.24 


9.1 x 10" 


13.9 


5.1 x 10" 


1.21 


26 


-88.5 


4.0 


20 


< 2.8 x 10 12 


1.8 x 10 12 


< 1.6 


<0.08 


1.5 x 10 12 


< 1.9 


4.5 x 10" 


< 0.31 


27 


-93.0 


4.0 


20 


9.9 x 10 12 


1.4 x 10 12 


6.9 


0.34 


2.5 x 10 12 


4.0 


6.0 x 10" 


0.80 


28 


-97.5 


4.0 


20 


1.4 x 10 13 


1.1 x 10 12 


12.6 


0.62 


2.0 x 10 12 


6.9 


6.0 x 10" 


1.15 


29 


-101.8 


4.0 


20 


1.1 x 10 13 


9.6 x 10" 


11.8 


0.58 


4.3 x 10 12 


2.6 


1.1 x 10 12 


0.48 


30 


-106.2 


4.0 


20 


1.1 x 10 13 


4.5 x 10" 


25.2 


1.24 


4.1 x 10 12 


2.8 


6.0 x 10" 


0.92 


31 


-111.6 


4.0 


20 


8.4 x 10 12 


< 1.6 x 10" 


>53.0 


>2.60 


8.6 x 10" 


9.8 


5.4 x 10" 


0.76 


32 


-116.0 


4.0 


20 


1.4 x 10 13 


< 1.6 x 10" 


>88.3 


>4.33 


4.3 x 10" 


32.8 


4.2 x 10" 


1.64 



Notes: All components have a rotation temperature of 2.7 K except components n c = 1 and 2 of C-C3H2 which have 15.0 and 3.0 K, respectively, 
and component n c = 2 of H 13 CO + which has 3.1 K. The emission components of the H 13 CO + model are the same as in Tableland are not listed 
here. 

8 Component number. * LSR velocity, the same for all four molecules except for 1 component of 13 CH + (2: 68.7 km s -1 ) , 2 components of c-C 3 H 2 
(1: 61.0 km s -1 , 2: 65.0 km s" 1 ) , and 1 component of H 13 CO + (2: 64.5 km s" 1 ). c Linewidth (FWHM), the same for all four molecules except for 
2 components of C-C3H2 (1: 15.0 km s _1 , 2: 11.0 km s -1 ) and for 1 component of H' 3 CO + ( 2: 12.0 km s -1 ). d Assumed 12 C/ 13 C isotopic ratio for 
HCO+ and CH + . ' Column density. f Column density ratio SH + / 13 CH + . * Column density ratio SH + / 12 CH + . * Column density ratio SH + /c-C 3 H 2 . ' 
Column density ratio SH+/HCO+. 



the other hand, this transition could partly contribute to the 
broad feature detected in emission between 58 and 125 km s _I 
(although the higher velocity half of that feature is uncertain due 
to the limited baseline quality). An additional model component 
in emission is indeed permitted by the 893 GHz spectrum and 
can even solve the disagreement at 62 km s _1 between the 
synthetic and observed spectra of the I0-O1 transition. If this 
emission component has a size smaller than ~ 20", then the T mt 
= 7cmb synthetic spectrum of the (-20.4", +1.4") CHAMP+ 
pixel at 683 GHz (shown in Fig. 13 will be unaffected and still 
provide a good fit to the observations. However, the fit to the 
spectrum of the central pixel at 683 GHz (see Fig. |5j will be 
worse, and a more elaborate radiative transfer modeling would 
be needed to find a consistent fit for this potential emission 
component. 



4.1.3. The 13 CH + absorption 

The 13 CH + 1-0 transition detected in absorption toward 
Sgr B2(M) was also modeled in the T lot = 7cmb approximation 
(see Fig. [2}:). We used the SH + model as a reference, keeping 
the velocity offsets and linewidths unchanged. The resulting 
fit, shown in red in this figure, is in good agreement with the 
observed spectrum, although it could certainly be somewhat 
better if the line widths were taken as free parameters for the 
fitting. 

The column densities resulting from this fit are listed in 
Table [3] as well as the SH + / 13 CH + and the SH + / 12 CH + column 
density ratio. The latter varies significantly from component to 
component, from 0.07 to 1.6, if we exclude component 32 in 
Table [3] at whose velocity the signal-to-noise ratio in the SH + 
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Table 4. Parameters of the best-fit LTE model of H 13 CO + and the LTE model of SH + fitted with the same velocity components as 
H 13 CO + . 







H 13 CO 


+ 




SH + 




n c 


1/ b 
Hsr 


A V 


7v 


Be 


ly 


p / 

"12 




(km s ') 


/i -1 \ 

(km s ') 


(cm 2 ) 




(cm -) 




i 1 \ 

W 


/"">\ 


(3) 


(4) 


(5) 


(o) 




1 
1 


/io n 
02. 1) 


1 o n 
12.0 


i oc w i nJ4 
3. /5 X 10 


on 
20 






2 


on n 


on n 
20.0 


c nn w i nl3 
5.UU X 10 


on 
20 






3 


104.0 


20.0 


/i nn w inl2 
4.00 X 10 


20 






4 


c 1 n 


1 o n 
12.0 


a nn v i nl3 
4.00 X 10 


on 
20 






c 

J 


/i/i c 
64.5 


1 o n 
12.0 


1 ^i/i w i nl4 
1.36 X 10 


on 
20 


i.y / x iu 


0.0/ 


c. 




q c n 
35. U 


o n 
/.U 


/i oo i nl 1 
4.82 X ID 


-1 n 
40 


< 2.o2 X 1U 


< U. 14 


-7 

/ 


n 

3U.U 


c n 
5.0 


H OQ w i nl 1 

/.83 X 10 


i n 

40 


^ o oo inl2 
< 2.o2 X 10 


< u.uy 




2D. 3 


1 n 
3,1) 


q n i w i nl 1 
3.U1 X 10 


i n 
40 


< 2.S2 X 10 


< 0.23 


y 


22.0 


c n 
5.0 


o w i nl2 
2.53 X 10 


/in 
60 


0.4j X 10 


0.05 


1U 


i o n 
i IX) 


c n 
5.0 


o nc w i nl2 
2.05 X 10 


/in 
60 


2.2j X 10 


n i o 
0. lo 


1 1 


1 1 n 
12.U 


/i n 
6.0 


1 o i w i nl2 
1.81 X ID 


/in 
60 


o oo -.^ 1 nl3 
2.o2 X 10 


n oc 
0.25 


i o 
12 


~7 /i 

/.o 


4.5 


1 on w i nl2 
1.21) X 10 


on 
20 


2.j4 X 10 


1.03 


13 


1 1 

1 . 1 


t-. o 
6.2 


/i 1 c v i nl2 
4.15 X 10 


on 
20 


l.OV X 10 


n on 
0.2U 


14 


-J.J 


"7 n 


^ no w i nl 1 
6.02 X 10 


on 
20 


1 n*7 \s 1 nl3 
Y.yi X 10 


l.ol 


15 


-13.5 


5.0 


i ni w inll 
3.1)1 X 10 


40 


1,55 X 10 


1.26 


lo 


1 o c 


c n 
5.0 


/;o w i nl 1 
6.62 X 10 


-1 n 
40 


o az \/ 1 nl2 
6.4j X 10 


n q 1 
0.3 1 


1 / 


O/i 0. 

-20.3 


c n 
5.0 


h oo w i nl 1 
1 .11 X 10 


A n 
40 


o a c \^ 1 nl2 
0.4j X 10 


n on 

o.2y 


1 O 


-31.0 


c n 
5.0 


o 1 o v i nl 1 
2.1 / X 10 


-1 n 
40 


o /i c \s i nl2 
6.4j X 10 


n n^ 


1 ("1 


-30. J 


-1 n 


1 on \y inll 
1.Z0 X 1U 


i n 
40 


o a z \y 1 nl2 
o.4j X 10 


1 TO 

1. /2 


20 


/tin 
-41.1) 


j.j 


1 co w i nl2 
1.5 / X ID 


-1 n 
40 


h c\a \s 1 nl2 
/.04 X 10 


n i 1 
0.1 1 


O 1 

21 


A Z O 

-45. / 


i n 
3.U 


c no w i nl 1 
6.U2 X ID 


-1 n 
40 


a oq \/ i nl2 
4.23 X 10 


n i o" 
0.1 / 


oo 
22 


/I O C 

-4y.j 


n 
J.U 


t fii inll 
j.Ul X 1U 


i n 
40 


a oq \y 1 nl2 
4.23 X 10 


n ^/i 
0.54 


O 1 

2 J 


c a n 
-54.0 


-1 n 
4.1) 


1 nn w i nl 1 
1.99 X 1U 


on 
20 


n ir\ \s 1 nl2 
V.30 X 10 


o on 

2.2y 


O 1 

24 


f(l O 

-5y.3 


-1 n 
4.1) 


c oo w i nl 1 
5. 11 X 1U 


on 
20 


^ o oo inl2 
< 2.o2 X 10 


< 0.24 


o c 
25 


-60. U 


c n 
J.U 


a o i \y inll 
4.21 X 1U 


on 
ZU 


, qo 1 nl2 
< 2.o2 X 10 


< 0.53 


0£ 

2o 


- /2.0 


c n 
5.U 


a o i w inll 
4.21 X 1U 


on 
20 


n o/i 1 nl2 
V.OO X 10 


1.15 


OO 

2 / 


O/i O. 
- /0.3 


-1 n 
4.1) 


A 01 v> 1 nl 1 
4.21 X ID 


on 
20 


o /i c \s 1 nl2 
6.4j X 10 


n no 


1 O 
25 


on i 
-S0.1 


-1 n 
4.1) 


o o i w inll 
2. /I X 1U 


on 
20 


"7 r\/i v/ 1 nl2 
/.04 X 10 


1 OO" 

1.2 / 


2 V 


O A O. 

-84.3 


c n 
5.U 


c i nl 1 
5. 12 X 1U 


on 
20 


1 o"7 w 1 nl3 
1.2/ X 10 


1 O 1 

1.21 




-88 5 


3 5 

j.j 


4.21 x 10" 




< 1 R? v 10 12 

^ Z..OZ, A. 1U 


<: W 


31 


-92.6 


3.5 


5.42 x 10" 


20 


8.45 x 10 12 


0.76 


32 


-97.2 


4.0 


6.62 x 10" 


20 


1.41 x 10 13 


1.04 


33 


-102.2 


5.0 


1.57 x 10 12 


20 


1.69 x 10 13 


0.53 


34 


-108.0 


4.0 


3.61 x 10" 


20 


1.13 x 10 13 


1.53 


35 


-112.0 


4.0 


4.82 x 10" 


20 


8.45 x 10 12 


0.86 


36 


-116.0 


4.0 


3.61 x 10" 


20 


1.41 x 10 13 


1.91 


37 


-120.0 


4.0 


1.20 x 10" 


20 


5.63 x 10 12 


2.29 


38 


-123.7 


4.0 


1.20 x 10" 


20 


5.63 x 10 12 


2.29 


39 


-130.4 


8.0 


1.99 x 10" 


20 


7.04 x 10 12 


1.74 



Notes: All components have a rotation temperature of 2.7 K except components n c = 1 to 4 of H 13 CO + which have 20.0 K and n c = 5 of H 13 CO + 
which has 3.1 K. " Component number. b LSR velocity, the same for both molecules except for one component of SH + (5: 64.0 km s~'). 
Linewidth (FWHM), the same for both molecules except for one component of SH + (5: 19.0 km s~'). d Column density. e Assumed isotopic ratio 
H 12 CO + /H 13 CO + . f Column density ratio SH+/HCO+. 



spectrum is very low and which is not detected in the 13 CH + 
absorption at all. 

The 13 CH + 1-0 spectrum shows an absorption feature at 
+ 172.6 km s . Although absorption in CO is detected at this 
velocity in our 3 mm line survey in the CO 1-0 and 13 CO 1-0 
lines (see Fig. [6j, this component is neither detected in H 13 CO + 
1-0, C-C3H2 2i2-loi, nor in SH + 1 1 — Oi (see Fig. [2J. Features at 
this velocity are usually attributed to the E xpanding Molecular 
Ring (EMR, e.g. iMorris & Serabvnlll996l) but +172.6 km s _I 
would refer to the far side of the ring, hence no absorption is 
expected and it is possibly more likely that this feature is an as 
yet unidentified (U) absorption line. 



We also list in Table [3] the parameters of an alternative 
H 13 CO + r, ot = Tqmb model fitted with the same velocity 
components as 0C3H2 (i.e. same linewidths and LSR 
velocities). This model is less good than the one that we 
described above and that we showed in Fig.|2^, but it can be used 
to roughly estimate the SH + /HCO + column density ratio (Col. 12 
of Table [3]). This ratio varies by about one order of magnitude, 
from 0.07 to 2.38. 
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Velocity for CO (km s ) 



-400 -200 




_C0(3-2) 
; SH+da-Oj) 



-200 







200 



400 



Velocity for SH (km s ) 



Fig. 6. a Spectrum obtained toward Sgr B2(M) at the frequency 
of the SH + Iq-Oi line with the APEX telescope in main beam 
brightness temperature. The red curve is the synthetic spectrum 
described in Sect. 14.11 the horizontal dotted line shows the 
continuum level, and the vertical dotted line marks the systemic 
velocity of the source (62 km s _1 ). The strong absorption lines 
are spiral arm components of CO 3-2. The lower axis gives the 
LSR velocity for SH + , while the upper axis refers to CO, with the 
same labeling as in panel b. b CO(1-0) spectrum obtained toward 
Sgr B2 as part of our line survey with the IRAM 30 m telescope. 
The horizontal dotted line indicates the zero level (after baseline 
removal because the observations were performed in position- 
switching mode with a far OFF position, which yields a very 
uncertain baseline level). In both panels, the dashed line marks 
the systemic velocity of the source in CO 3-2 and the dot- 
dashed line indicates the CO velocity +172.6 km s at which 
an extra absorption component appears in the 13 CH + spectrum 
(see discussion in 34.1.31 



4.1 .4. Modeling the HCI absorption only from the Sgr B2 
region 

The HCI and H 37 C1 transitions detected in absorption toward 
Sgr B2(M) were also modeled in the r iot = 7cmb approximation 
(see Fig. [8]). We used the SH + model as a reference, keeping 
the velocity offsets and linewidths unchanged, except for the 
main absorption component associated with Sgr B2(M) itself. 
The resulting fit is shown in red in Fig. [8] and the fit parameters 
are listed in Table [5] Only the main absorption component is 
clearly detected for both isotopologues. For all other velocity 
components but one, we used a column density of 2.8 x 10 12 
cm -2 , represented by the red line in Fig. [HJ to get a rough 
estimate of the column density upper limits. Only one additional 
velocity component (n c — 6; ~ 17 km s -1 ) marginally shows 
absorption in the main isotopologue HCI at 2.8 times this value. 
Due to the compact hyperfine structure, the line width of the 
main absorption component (8 km s ) is somewhat uncertain, 
but our modeling suggests that it is significantly smaller than for 
the other species. Finally, for Sgr B2(M) we derive an isotopic 
[H 35 C1/H 37 C1] ratio of ~ 4, consistent with the terrestrial value. 



Rest frequency (GHz) 
893.5 893 



35 



30 



25 = 



20 - 



" SH + (2 1 - 


1,) 


- 










-200 



-100 







100 



200 



Velocity (km s ) 



Fig. 7. Spectrum obtained toward Sgr B2(M) at the frequency 
of SH + 2 1 - 1 1 with the APEX telescope in main beam brightness 
temperature. The red curve is the synthetic spectrum described 
in Sect. 14. II i.e., showing no detectable emission, the horizontal 
dotted line shows the continuum level, and the vertical dotted 
line marks the systemic velocity of the source (62 km s _1 ). The 

absorption line appearing at 105 km s _I is the Jk„k~ = 1 11 ~ 

Oon g round-state line of HDO at 893638.7 MHz dComito et al.1 

Rest frequency (GHz) 
626.5 626 625.5 625 



^^^^^ 


| , , , , : _ 




h 37 ci(i-o) : 

, i , i l , i , i l , i I i l , i r 


-200 -100 





100 


-200 -100 





100 



Velocity for HCI (km s ] 



Velocity for H CI (km s ) 



Fig. 8. Spectrum obtained toward Sgr B2(M) around 625 GHz 
with the APEX telescope in main beam brightness temperature. 
The red curve is the synthetic spectrum described in Sect. 14.11 
the horizontal dotted line shows the continuum level, and the 
left and right vertical dotted lines mark the systemic velocity 
of the source (62 km s _1 ) for HCI and H 37 C1, respectively. The 
strong emission close the HCI is likely due to a SO2 line at 
626087.3 MHz. 

Table 5. Parameters of the best-fit LTE model of HCI and H 37 C1. 



Molecule 


K 




AV C 


N d 






(kms- 1 ) 


(km S" 1 ) 


(cm- 2 ) 


(1) 


(2) 


(3) 


(4) 


(4) 


HCI 


2 


65.0 


8.0 


1.90 x 10 14 




6 


17.0 


5.0 


7.60 x 10 12 


H 37 C1 


2 


65.0 


8.0 


4.75 x 10 13 



Notes: " Same component numbering as in Table [3] b LSR velocity. 
Linewidth. d Column density. 



4.2. Summary of modeling 

To visualize the results of our modeling presented in Tables [3] 
and|U we show, in Fig.|U the column densities and abundances 
that we derived for the species considered functions of LSR 
velocity. 
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5. Discussion 

5.1. Column densities and abundances 

5.1.1. Caveats 

We would like to precede our discussion of the SH + , 13 CH + , 
and HQ abundances in the spiral arm clouds with some 
cautionary notes. First, the assignment of velocity features to 
Galactocentri c distance is by no means clear-cut. For example, 
ISofuel (l2006l) . from a careful analysis of the longitude-velocity 
diagram of the GC region, assigns certain velocity features to 
peculiar components of the complex Galactic center velocity 
field, namely the GC molecular ring and the GC expanding 
ring, both within a Galactocentric radius, R, of 600 pc of the 
GC. In addition to several spiral arms, he invokes expanding 
rings at R = 3 and 4 kpc. Some of Sofue's conclusions are at 
variance from the LSR velocity-location assignment we present 
in 34. II Second, SH + and 13 CH + may in principle be associated 
with predominantly atomic or molecular diffuse clouds. For 
meaningful abundance calculations, it is thus necessary to have 
values for both the atomic hydrogen column density, A^(HI), 
and the molecular hydrogen column density, Af(H 2 ). Actually, 
the clouds from which we see absorption may have complex 
structure, e.g., a Giant Molecular Cloud region and a diffuse 
atomic envelope (see 35.2b and/or we may be dealing with more 
than one cloud along each line of sight. Certainly, the HI spin 
temperatures are in the range of values found for diffuse clouds 
and significantly higher than the canonical GMC values of 20- 
30 K. 

Two final, more technical caveats relate to the way we 
derive column densities. For the 3 mm Sgr B2(M) spectra 
for H 13 CO + and c-C 3 H 2 we have constructed full XCLASS 
models considering all the absorption from these molecules plus 
emission from them and all other molecules. In contrast, for 
the submillimeter spectra, we have not done this. While even 
the 3 mm spectra may be contaminated by emission features 
inadequately treated by the model, for the submm absorption 
spectra this is a much greater source of uncertainty for line 
assignments and intensity predictions, making, strictly speaking 
all the column densities we determine lower limits. We do at 
present not feel confident to construct a comprehensive model 
for the submm spectra. Therefore, we have just marked the 
positions of stronger lines in the spectra by the species from 
which the originate. Similarly, if, as discussed in ^4] excitation 
by the ambient radiation field is non-negligible, then T rot > 
Tcmb, which would make our column densities underestimates. 

Finally, our XCLASS-based method makes it difficult to 
derive statistical errors and, thus, upper limits rigorously in a 
X 2 sense. Rather, upper limits are determined by assuming a 
column density, calculating a model spectrum and then "eye- 
balling" whether it is consistent with the data. 

5.1.2. Atomic and molecular hydrogen column densities 

The column density of atomic hydrogen, A^(HI), can be 
calculated from measurements of the 21 cm HI line's optical 
depth, t(HI), at t he intervening clouds' velocities (e.g., 
iGarwood & Dickevl 1 19891) and a determination of their spin 
temperatures, Ts, as described bv lLazare 1 (119751) . The method 
invokes a combination of interferometric observations to, 
first, directly determine the optical depths of the absorption 
components and, second, single dish observations to measure 
the (optically thick) line's brightness temperature. As described 
in a forthcoming paper, we have performed the necessary single 



dish observations with the Effelsberg 100 meter telescope and 
obtained an optical depth spectrum (t(HI) vs. V/ Jr ) measured 
toward Sgr B2(M) with the NRAO Very Large Array from 
C. Lang (pers. comm.). We note o ur values for jV(H I are 
generally lower than those derived bv lVastel et alj (Hp02). This 
is partially due to their assumption of T s — 150 K for all velocity 
components. 

A determination of Af(H 2 ) is more difficult, as selection 
of a "universal" trace molecule for molecular hydrogen is not 
straightforward. For meaningful modeling, such a molecule's 
lines would be required to be subthermally excited, meaning 
that their relative level populations be determined by Tcmb - This 
means, for example that, due to their small dipole moments 
and resulting low critical densities, isotopologues of CO are not 
suited. A possibl e choice for an H? trac er, which we adopt here, 
might be HCO + . iLucas & Lisztl (Q996) from their study of this 
molecule's 89 GHz 7=1-0 ground-state transition, conclude 
its abundance (relative to H 2 ) to be "~ 3-6 x 10~ 9 across a very 
broad range of extinction", meaning between hydrogen column 
densities ~ 0.1 and 4 x 10 20 crrT 2 . These authors base this 
finding, first, on the very tight correlation between the HCO + 
and the OH column density and, second, on the fact that OH has 
the same abundance (10~ 7 ) in diffuse as in dense dark clouds. 
While their own absorption data only sample hydrogen column 
densities between 10 20 and 2 x 10 21 cm, -1 , they find that the 
correlation holds at 200 times higher values. Here we adopt 
5 x 10~ 9 for the HCO + /H 2 abunda nce ratio. This value is of 
the same magnitude as that listed bv lvan Dishoeck et alj d 19931) 
for (high column density) regions of low and high mass star 
formation. 

Since HCO + absorption toward Sgr B2 is optically thick even 
for the non-Sgr B2 region velocity features, we use, as described 
in 34.11 H 13 CO + spectra and obtain the main-isotope's column 
density by scaling with the appropriate 12 C/ 13 C ratios given in 
column 4 of Table [3] 

Our estimates for jV(HI) and A^(H 2 ) for the absorption 
components are listed in Table [6] For simplicity's sake and 
motivated by our optical depth fitting of the HI data and also 
to allow comparison with other studies, we have grouped the 
velocity components into seven intervals. The HI and H 2 column 
densities and HI spin temperatures were derived as described 
above from HI and H 13 CO + data fixed to values compatible 
with the values given in the table. The column densities of SH + 
and CH + , presented in Table [7] were obtained by summing over 
the individual velocity components in Table [3] and multiplying 
the CH + column densities by the 12 C/ 13 C ratios given for each 
component in that table. 

5.2. The origin of the 13 CH + and SH + absorption 

Since it is a priori not clear whether the observed SH + , CH + 
and HC1 molecules are only associated with either atomic or 
molecular gas or both, when calculating abundance ratios, we 
might want to refer to the total column density of hydrogen 
nuclei, N(H) = N(Hl) + 2 N(H 2 ). Inspecting Table we find 
that over all velocity ranges showing absorption, A^(H 2 ) greatly 
exceeds A^(HI), by factors between 13 and 133, approaching 
the value for Sgr B2 (138). This means that the absorbing 
medium is dominat ed by GMC gas. Thi s is reinforced by the 
densities derived by iGreaves et alJ (1 19921) from LVG modeling 
of the two lowest rotational lines of CS, which, between 0.4 and 
2 x 10 4 cirT 2 , are actually even higher than values in most of a 
GMC's volume. 
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Table 6. H i and H2 column densities of the absorption components. 



T 

La Del 


Velocity range 


\r IT b 


-r c 
1 s 




Al e 


Abs. comp. 




(km s _1 ) 


(cm -2 K" 1 ) 


I'K ) 


(cm - ) 






(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


A 


95.3-38.6 


9.0(1.8) x 10 19 


43(5) 


3.9(0.9) x 10 21 


5.4 x 10 23 


[2] 


B 


23.2-12.9 


1.5(0.5) x 10 19 


65(10) 


1.0(0.4) x 10 21 


5.5 x 10 22 


[5,6] 


C 


-12.9 --33.5 


8.5(2.0) x 10 18 


93(10) 


7.9(2.0) x 10 20 


1.4 x 10 22 


[12, 13, 14, 15] 


D 


-36.1 --5 1.5 


1.6(0.4) x 10 19 


20( 3) 


3.1(0.9) x 10 20 


2.0 x 10 22 


[17, 18, 19] 


E 


-54.I --64.4 


3.0(1.0) x 10 18 


55( 6) 


1.7(0.6) x 10 20 


2.6 x 10 21 


[21] 


F 


-69.6- -85.0 


2.6(0.7) x 10 18 


80( 8) 


2.1(0.6) x 10 20 


7.1 x 10 21 


[22, 23, 24, 25] 


G 


-90.2- -123.7 


2.8(0.7) x 10 18 


43(7) 


1.2(0.4) x 10 20 


1.6 x 10 22 


[27, 28, 29, 30,31,32] 



Notes: the numbers in parentheses are the uncertainties. " Velocity range over which the Hi column density is computed. * Ratio of Hi column 
density to spin temperature, obtained from the VLA data. c H 1 spin temperature derived from the Effelsberg data. d H 1 column density. e Hi column 
density derived from the H 13 CO + column densities assuming the isotopic ratios listed in Col. 4 of Tableland an HCO + abundance of 5 x 10~ 9 
relative to H2. s List of the absorption components used to compute the H2 column densities, with the same numbering as in Col. 1 of Table[3] 



Both the total hydrogen nuclei column density, N(H), and the 
H2/HI ratio are muc h larger than the val ues found in translucent 
and diffuse clouds. iGredel et al.1 (1 19921) in their study of the 
archetypical translucent cloud toward HD 210121 find N(H2) = 
8xl0 20 ctrT 2 andAf(HI) = l-5xl0 20 cirT 2 for the molecular and 
atomic hydrogen absorption column densities, respectively. In 
comparison, the classic diffuse cloud line of sight toward g Oph 
has comparable HI and H2 column densities of * 5xl0 20 cirT 2 or 
a column density of hydrogen nuclei of A^(H) = 1.5 x 10 21 cm -2 
(I van Dishoeck & Blackll 19861) . 

As discussed above, the absorbing material is dominated 
by GMCs. There is, however, nothing that does imply that the 
observed CH + and SH + is associated with gas within a GMC, 
although little is known about the abundance of both species 
in such objects. Model calculations predict significant column 
densities in lower-Ay portions of dense (> 10 5 cm s~ 3 ) PDRs 
illum inated by intense ultraviolet fields dSternberg & Dal garno 
Il995l) . Even more extreme conditions (XDRs) are required to 
produ ce significant amounts of SH + in dense regions ( Abe l et al.l 
2008). However, both dense P and XDRs have small volume 
filling factors, making them unlikely absorbers. These findings 
point to an (at least) two component picture for most of the 
absorption components, i.e., a diffuse envelope traced by HI, 
CH + and SH + abs orption plus a dens er core region, a view that 
is also invoked by Vast el et al.1 (|2002|) to explain the absorption 
they observed with ISO in the C + and O fine structure 
transitions. 

In Fig. [9] we plot CH + column density vs. hydrogen column 
density, A^(HI) as defined above, for the seven velocity ranges 
in Table [6] (coded A-G). For each CH + column density entry 
two data points are plotted, both taken from the table. The left 
one (open squares) assumes that the CH + is only associated with 
atomic hydrogen, where the right one (filled squares) assumes 
it is only associated with molecular hydrogen or, essentially, all 
the hydrogen since A^(H2) » A^(HI) for all velocities. We also 
plot (as circles) data for the two classic diffuse clouds against £ 
Oph and f Per. For these, 7V(H 2 ) ~ iV(HI) ~ 5 x 10 20 cm 2 . 
In addition, we present data for the translucent cloud toward 
HD 210121 where, again, two points assuming that the CH + 
is only associated with atomic (left) or molecular gas (right) 
are plotted. All values for the latt er three sources, jV(H2), 
WHD, and Af(CH + ), are taken from Ide Vries & van Dishoeck! 
( 1988). The circles represent the absorption data Gredell dl997) 
determined toward a sample of southern OB associations. For 



these lines of sight, which have higher extinctions, between 0.36 
and 4.2, (mostly) corresponding to translucent clouds, we cannot 
discriminate between atomic and molecular hydrogen columns 
as Gredel lists A^(CH) + vs. Ay, which contains contributions 
from both. To convert his Ay values to hydr ogen column 
densit ies we use the results of the classic study of iBohlin et al.l 
(1978), who found for low values of visual extinction (0.02 < 
Ay < 0.6) a linear relation bet ween Ay, and hy drogen column 
density, which as first argued by Dickmanl (1 19781) . also holds for 
significantly higher Ay values. It is con sistent with a relatio n 
derived from modeling of x-ray data (iGliver &Ozell |2009): 
JV(H(cm 2 )) * 2 x 10 21 Ay(mag). 

5.2.1 . Comparison with individual diffuse clouds 

While we advocate above a diffuse cloud origin of the CH + 
absorption, we note that, apart from one velocity range (B, 12.9— 
23.2 km s _1 ), our CH + column densities are much higher than 
found in diffuse clouds and mostly higher than in the translucent 
cloud sample. However, the velocity ranges over which the 
values were determined, are mostly much larger than the line 
widths found in diffuse or translucent clouds (and even GMCs). 
For this reason, we plot in the lower panel of Fig. |9]the average 
column densities that a cloud with "typical" line width AV typ 
would have. This means we have divided for each of the velocity 
intervals A-G, both, CH + and hydrogen column densities by the 
ratio of the width of the velocity interval and AV typ . We chose 
AV typ to be 4.3 km s FWHM, whi ch is the mean of t he sample 
of diffuse lines of sight for which ICrane et al ] (119951) present 
high spectral res olution CH + d ata0 This value is consistent with 
that discussed bv lGredeUd 19971) for his translucent cloud sample. 
The thus calculated values correspond to a "typical" cloud with 
averaged properties. 

Inspecting the CH + column densities for these typical 
clouds we find values very similar to those found in diffuse 
and translucent clouds. We note that the ratio of molecular 
to atomic hydrogen column densities is in most cases higher 
than that found in translucent and much higher than in diffuse 
molecular clouds. The reason for this is almost certainly sight 
line superposition of the CH + -containing diffuse/translucent and 
denser molecular clouds, as discussed above. 

8 The Doppler parameter, b, listed by this and numerous other optical 
absorption studies and FWHM line width are related via AV 1/2 = 
l^hUb a 1.665 b. 
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Fig. 9. CH + column density versus hydrogen column density. 
The ordinate gives the CH + column density for the seven 
different velocity intervals listed in Table [7] The capital letters 
designate the entries in table (top to bottom). The open circles 
re present data tak en from the diffuse/translucent cloud sample 
of iGredell (fl997h . The squares connected by lines mark the 
CH + column densities determined in this work. For the open 
squares it is assumed that all the absorption is from atomic gas; 
the abscissa corresponds to the HI column density (Table [6]), 
whereas for the filled squares it is assumed that all the absorption 
is from molecular gas; the abscissa corresponds to the H2 column 
density (Table [6]). For the upper panel, column densities were 
calculated from the absorption over the whole velocity range 
of each interval, whereas in the lower panel average column 
densities normalized to a typical line width of 4.3 km s~' were 
used (see text). 



5.3. Comparison - SH + 
CH+/SH+ ratio 



and other S-bearing species and the 



The values we derive for the SH + column density, A^(SH + ), for 
the non-Sgr B2 region velocity ranges go up to 7 x 10 13 cm -2 
(see Table |7). Considering again the differences in line widths, 
these values are comparable with the upper limits found for 
optical SH + absorption from diffuse and translucent clouds ( i.e. 
1 x 10 13 cm" 2 ; see lL2b . 

The column densities of other sulfur bearing species in 
the velocity r anges of Table [7] have been determined by 
IGreaves et all (1 19921) . who studied CS, and iTieftrunk etall 
dl994j) (SO, H 2 S, and S0 2 ). Comparing with our values for SH + 
(Table [6), we find an [SH + /CS] abundance ratios between 0.5 
and 3, [SH + /SO] between 1.6 and > 55, and [SH + /S0 2 ] > 2. 
For the [SH + /H 2 S] ratio we find values > 2, except for the -33.2 
to -50.1 km s _1 velocity range, for which we find the low ratio 



of 0.3. We note that the velocity of this absorption places it in 
the Norma spiral arm. From Table [6] we see that this range has 
the highest H 2 column density and the highest H 2 /HI ratio of all 
the (non-Sgr B2 region) absorption components and resembles 
most closely a GMC or GMC core. This reinforces SH + 's nature 
as a tracer of more diffuse gas whose abundance drops in denser 
environments. 

As discussed by ITieftrunk et "all (11994 . SH + and H 2 S are 
chemically related, as H3S + , from which the latter forms via 
dissociative recombination is formed by the reaction 

SH + + H 2 -> H 3 S + + hv. 



With respect to chem ical models dMillar et alj 119861: 
iPineau des Forets et al.1 fl986). it is also worth noting that, as 
can be seen from Table □ the SH+/CH+ ratio, [SH+/CH+], 
which is 0.56 for Sgr B2(M) only varies between 0.11 and 

0. 26 for all but one of the non-Sgr B2 absorption features. 
These values are consistent with the upper limits on this ratio 
(0.12-2) for diffuse and translucent lines of sight (see jsll.2l i. 
The component with the higher ratio, * 0.9 is in the -90.2 
to -123.7 km s" 1 range, which correspond to peculiar velocity 
gas associated with the GC. Interestingly, this velocity range 
has also by far the highest SiO colu mn density measure d 
for the various velocity components by IGreaves et alj {1996). 
These two findings may have a common (shock-wave) origin 
as an incre ased SiO abundanc e is commonly ascribed to shock 
chemistry dSchilke et alj|1997l) . 

5.4. Limits on HCI absorption in the diffuse material 

The HCI column density, TY(HCl), we derive for the gas 
associated with Sgr B2(M) (see Table |5]l is i n good agreement, 

1. e. 16 % higher, with the value derived by IZmuidzinas et alJ 
(fl995l) . 

For the velocity components corresponding to the diffuse 
lines of sight we derive an upper limit on the HCI column 
den sity, jV (HCl), of 8.4 x 10 12 cm~ 2 (see Table and Fig. 
in 34.1.41 1. Since 7Y(H) ranges from 0.5-2.4 x 10 22 cnT 2 we 
calculate an HCI to hydrogen nuclei abundance ratio [HC1/H] 
of 0.6 - 3 x 10 -10 . Our upper limit is much c oarser than the 
(2.7 ± 1.0) x 10 11 cnr 2 iFederman et all (1 19951) determine for 
Af(HCl) from an UV absorption line toward f Oph. 

6. Conclusions 

We have reported submillimeter wavelength absorption from 
rotational ground-state transitions of 13 CH + , H 35 C1, H 37 C1, and, 
for the first time, of the SH + molecule toward the strong Galactic 
center continuum source Sgr B2(M). From I3 CH + and SH + we 
detect absorption at the systemic LSR velocity of Sgr B2(M) 
and in addition over a wide range of lower velocities belonging 
to intervening material. For the intervening clouds we have 
derived molecular hydrogen column densities from optically thin 
H I3 C0 13 emission and atomic column densities from HI 21 cm 
line data. For all the clouds we find that the former is larger 
than the latter by factors between 15 and 130. If we "scale" 
the observed velocity ranges of the individual clouds to the line 
widths of single diffuse and translucent molecular clouds in the 
Solar neighborhood, we find similar 13 CH + column densities 
in both types of clouds. Most likely the reason for this is that 
the observed 13 CH + is indeed mostly associated with atomic 
material surrounding molecular material. The [SH + /CH + ] ratio 
is consistent with upper limits from optical data. Our upper 
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limits on the HQ column density in the absorption components 
not associated with Sgr B2 are consistent with but coarser than 
measurements of UV absorption from this molecule. 



Table 7. Column densities of the absorption components and abundances relative to H. 





SH 


+ 


CH 


+ 




SO 


H 2 S 


so 2 




CS 


Velocity range 


N 


X H 


N 


X H 


N 




N 


X H 


N 




N 




(km s -1 ) 


(cm ) 




(cm - ) 




(cm ) 




(cm- 2 ) 




(cm- 2 ) 




(cm ) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


95.3-38.6 


2.0(14) 


1.8(-10) 


3.6(14) 


3.3(-10) 


















23.2-12.9 


3.2(13) 


2.9(-10) 


5.7(13) 


5.2(-10) 


















-12.9 --33.5 


3.4(13) 


1.2( -9) 


3.1(14) 


l.l(-8) 


3.6(12) 


1.0(-10) 


1.5(13) 


5.2(-10) 


< 2.1(13) 


<7.3(-10) 


1.6(13) 


5.5(-10) 


-36.1 - -51.5 


1.5(13) 


3.9(-10) 


9.2(13) 


2.3( -9) 


9.4(12) 


2.0(-10) 


4.8(13) 


1.2( -9) 


< 1.1(13) 


<2.7(-10) 


2.8(13) 


7.1(-10) 


-54.1 --64.4 






1.4(13) 


2.6( -9) 


















-69.6- -85.0 


3.9(13) 


2.7( -9) 


1.5(14) 


1.0( -8) 


< 7.0(11) 


<5.0(-ll) 


< 4.2(12) 


<2.9(-10) 


< 3.7(12) 


<2.6(-10) 


1.3(13) 


9.1(-10) 


-90.2- -123.7 


6.9(13) 


2.2( -9) 


7.9(13) 


2.5( -9) 


5.1(12) 


2.0(-10) 


< 4.5 (12) 


< 1.4(-10) 


1.9(13) 


6.1(-10) 


4.4(13) 


1.4( -9) 



Notes: A' is the column density of molecule Y and X H the column density ratio [Y]/([Hi]+2x[H2]). The column densities of SO, HtS, and SO2 are from Tieftrunk et al. (1994) and those of CS 
are from Greaves et al. (1992). The velocity ranges are the same as in Table [6] The list of absorption components used to compute the SH + and CH + column densities is given in Table [6] The 
components for which the column density is an upper limit in Tablef3]were ignored. x(y) means x x 10 v . 
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